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Abstract 

We consider solid surface scattering of molecules that were subject to strong non-resonant ul- 
trashort laser pulses just before hitting the surface. The pulses modify the rotational states of the 
molecules, causing their field free alignment, or a rotation with a preferred sense. We show that 
field-free laser induced molecular alignment leads to correlations between the scattering angle and 
the sense of rotation of the scattered molecules. Moreover, by controlling the sense of laser in- 
duced unidirectional molecular rotation, one may affect the scattering angle of the molecules. This 
Q-t provides a new means for separation of mixtures of molecules (such as isotopes and nuclear-spin 

isomers) by laser controlled surface scattering. 

PACS numbers: 
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I. INTRODUCTION 



Laser control of molecular rotation, alignment and orientation has received significant 
attention in recent years (see e.g. reviews jlj). Interest in the field has increased, mainly 
due to the improved capabilities to manipulate the characteristics of the laser pulses (such 
as time duration and temporal shape), which in turn lead to potential applications offered 
by controlling the angular distribution of molecules. Since the typical rotational time scale 
of small molecules is 'long' (~10ps) compared to the typical short pulse duration (~50fs), 
effective rotational control and manipulation are in reach. During the last decade, coherent 
rotational dynamics of pulse-excited molecules was studied 
giving rise to enhanced alignment were suggested 
manipulations, such as optical molecular centrifuge [5| and alignment-dependent strong field 
ionization of molecules [6fl, were demonstrated. Selective rotational excitation in bimolecular 



2] , and multiple pulse sequences 



and realized experimentally [4|]. Further 



mixtures was suggested and demonstrated in the mixtures of molecular isotopes [7( and 
molecular nuclear-spin isomers [8J. Aligned molecules were also shown to be useful for 



other applications, including: generation and modification of ultrashort 
control of high harmonic generation [lo| . amplification of terahertz pulses 
of molecular structure 



n 



ight pulses 



ll|, elucidation 



121 ]. rotational spectroscopy [13], manipulation of chemical reactions 



14j, and others. Finally, it was demonstrated that these new methods for manipulation of 



molecular rotation can also be used to modify the motion of molecules in inhomogeneous 



fields, such as focused laser beams 
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16] , or static electric [17] and magnetic [18] fields. 



Modification of the molecule-surface scattering and molecule-surface reactions by exter- 
nal fields of different nature is a long-standing research problem. In particular, the effect 
of molecular orientation by a static electric hexapole field on the scattering process was 
investigated in detail []j]]. Laser control of the gas-surface scattering process was achieved 
using multiphoton ionization of the impinging molecules by long laser pulses of variable po- 
larization 20[], and the possibility of controlling molecular adsorption on solid surfaces using 
ultrashort laser pulses was discussed [21]. 

In a recent paper [22], we investigated two possible schemes for modifying molecular 
scattering from flat surfaces. In the first scheme, we considered excitation of unidirectional 
molecular rotation via the combined action of a short laser pulse on molecules in a molecular 
beam, followed by a surface collision of these molecules. In the second scheme, we suggested 
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exciting a molecular beam, consisting of two molecular species, by two properly delayed 
laser pulses, in order to provide selective rotational excitation of one of the species. We have 
shown that these two species are then scattered differently from a solid surface, in a way 
that the scattered beam is enriched in one of the species, for specific scattering angles. 

The current paper extends the model of molecule-surface scattering that we used in 22]. 
We treat here surfaces that are not flat, but corrugated on different length scales. One limit 
includes surfaces whose corrugation period is larger than the size of the molecule. The model 
we developed according to this limit is based on the model in 23). This model is referred 
to throughout the paper as the corrugated surface model. This limit corresponds to highly 
corrugated surfaces, such as Pt(211) [23j, or LiF(OOl) |24|. Another limit includes surfaces 
with corrugation period smaller than the size of the molecule, or relatively flat surfaces. Here 
the model is based on the work of [25]. This model is referred to below as the surface with 
friction model. This limit corresponds to relatively flat surfaces, such as Ag(lll), where 
mechanisms other than surface corrugation exist for producing forces tangent to the surface, 
such as: tangentially directed phonons, creation of electron-hole pairs, etc. 25]. 

In this paper, we continue to investigate the ability to modify and to control the molecule- 
surface scattering processes by using ultrashort laser pulses. These laser pulses (mainly 
femtosecond laser pulses) align the molecules or set them in some specific rotational motion, 
before the molecules hit the surface. The paper is organized as follows: In Sec. [Til we present 
the model of molecular scattering from a corrugated surface. The classical model of interac- 
tion between an ultrashort laser pulse and a linear rigid molecule is briefly summarized in 
Sec. IIII1 Then, in Sec. IIVI it is shown that a laser pulse aligning the molecules before they 
scatter from a corrugated surface, leads to a correlation between the direction of scattering 
of the molecules and their sense of rotation. This correlation also depends on the direc- 
tion of polarization of the laser. The two-dimensional and the three-dimensional models of 
molecular scattering from a surface with friction is next presented in Sec. [V] Accordingly, in 
Sec. IVIl the scattering of molecules from a surface with friction, set in rotation in a specific 
sense by laser pulses, is investigated. It is shown, that a correlation exists between the sense 
of rotation of the molecules impinging on the surface, and their scattering angles. Finally, 
in Sec. IVIII we conclude and summarize. 
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II. A MODEL OF MOLECULAR SCATTERING FROM A CORRUGATED SUR- 
FACE 



Here we develop a model and analyze the scattering of a molecule from a surface that is 
corrugated at a length scale larger than the size of the molecule. This model can qualitatively 
describe the scattering process from surfaces, such as Pt(211) 

The model developed here is an extension of our model 



23. 

22j from flat to corrugated 



surfaces. Our treatment of corrugations is based on the model of Tully [23] . The diatomic 
molecule is treated as a rigid dumbbell {26 1. This dumbbell collides with a corrugated, 
but frictionless, surface, which is composed of hard cubes on identical springs, see Fig. [TJ 



Each cube represents one (or several) of the surface atoms [23j, |24|, We assume that 



the corrugation is one dimensional, and that its period is larger than the molecular bond 
length, such that the molecule collides with a single cube that is oriented in a definite 
angle. We also assume that the cube oscillates with an amplitude determined by the surface 
temperature. This hard cube model provides a simple way of adding surface phonons to 
the process. For the sake of simplicity, we also assume that the cube is much heavier than 
the molecule, so that its velocity does not change as a result of the collision. By moving to 
the reference frame attached to the cube, one reduces the problem to the molecular collision 
with a motionless hard wall. In the moving coordinate system, the molecular total energy 
(translational+rotational) is conserved, but it can be redistributed between these two parts 
as a result of the collision. 

What about the translational linear momentum of the center of mass of the molecule in 
the moving coordinate system? The momentum component that is locally perpendicular to 
the surface is not conserved, because the vibrating cube exerts forces on the molecule in 
this direction during the collision. On the other hand, there are no forces applied in the 
direction that is locally parallel to the frictionless surface, and therefore the linear momentum 
locally parallel to the surface is conserved. It follows, that because of the corrugation, the 
momentum that is globally parallel to the surface is not conserved, contrary to the case of a 
flat frictionless surface. 

Using energy and angular momentum conservation laws (angular momentum is conserved 
only in the coordinate frame centered at the position of the colliding atom), we find analytic 
expressions for the translational and the rotational velocities of the dumbbell molecule after 
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the collision. These velocities depend on the local velocities before the collision, and on 
the local angle between the dumbbell and the face of the cube at the moment of collision. 
Finally, we transform the velocities back to the laboratory coordinate frame. 

It is important to emphasize that our model treats the molecule as a dumbbell, although 



it is also possible to treat it as an ellipsoid 



24|. In our model, differently from the model 



used in [24|], we include three-dimensional rotation and multiple collisions of the molecule 
with the surface. 

In the next subsections, we first define the coordinates and the surface corrugation in 
subsection III Al Then, we consider the effects of the surface-cube vibration on the collision 
in subsection III Bl In subsection III C\ we provide the results for the molecular velocities 
after the collision as a function of the velocities before the collision and of the molecular 
orientation. 



A. The coordinates and the one-dimensional surface corrugation 

We treat a homonuclear diatomic molecule as a massless stick of length r e (equal to the 
bond length), with two balls, each of mass m, attached to its ends. To describe the motion 
of the molecular center of mass, we define the Z'- and X' -R.XGS clS globally perpendicular and 
parallel to the surface, respectively. We also define the Z- and X locally perpendic- 

ular and parallel to the surface, respectively, see Fig. [TJ The local molecular orientation is 
given by the vector r = 0.5r e (sin#cos<^, sin#sin<^, cos#). Notice that this molecular orien- 
tation depends on the orientation of the cube on which the molecule falls. Here the polar 
angle between the molecular axis and the Z-axis is 6, ranging from to n, and ip is the 
azimuthal angle, with values from to 2ir. The global translational velocity components of 
the center of mass are denoted by V' z and V' x . The first component is defined as positive 
if the molecule moves upwards, and the second one is positive if the molecule moves to the 
right. The local translational velocity components (in the unprimed coordinate system) are 
denoted by Vz and Vx- The linear rotational velocity is denoted by v, and is equal to dv/dt. 

The incident and the scattering angles are measured with respect to the Z'-axis, and are 
positive in the clockwise direction and negative in the counter-clockwise direction. The angle 
of incidence of the molecular beam is denoted by 8- mc , and is always chosen as non-positive. 
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FIG. 1: Schematic drawing of the dumbbell - hardcube corrugated surface model. The surface is 
a two-dimensional sinusoidal sheet composed of small hard cubes of mass M. The period of the 
sine corrugation is c and its amplitude is a. The cubes are held by springs that have a spring 
constant of Muj 2 . They oscillate in a direction perpendicular to the local surface tangent, denoted 
by the Z-axis. The inclination of the Z axis from the global surface normal (shown by the dashed 
line) is denoted by a. The molecule is a massless stick with two balls of mass m on its ends. The 
molecular orientation of one of the atoms is given by the vector r, and the linear velocity of this 
atom is v. The incident molecule approaches the surface at an angle of incidence of #i nc , which is 
always negative. The global coordinate system X'Z 1 is shown on the right. 

The angle between the local normal to the surface Z and the Z'-axis is denoted by a. The 
velocities in the primed coordinates can be expressed as a function of the velocities in the 
unprimed coordinates and the angle a: 

V x \ I cos a — sin a \ [V'x 
Vz J \ sin a cos a / \ V z 

Following the work of Tully [23] , we choose the simplest sinusoidal one- dimensional cor- 
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rugation of the form: 

Z'(X') = acosf— ) , (2) 



V c 

where a is the amplitude of the corrugation, and c is its period. Scattering out of the X' Z'- 
plane is not possible in the framework of this model, due to the absence of corrugation along 
the y'-direction, and to the absence of a local friction. This treatment corresponds to the 



281 ] . The model used in Sec. [V] incorporates out-of-plane 



measurement of in-plane scattering 
scattering due to frictional forces parallel to the surface. 

The tangent of the angle of the local surface normal a is given by: 

dZ' 2na (2itX'\ 
tana = -— = — su^—j . (3) 

The maximal local normal angle is a max = arctan (27ra/c), which is referred to as corrugation 
strength. A molecule incident on the surface at an angle ^inc, hitting a cube inclined at a, has 
an effective incident angle of #i nc — a. We will always assume that no "shadowing" occurs, 
which means that 9 mc — a max > —ti/2. We also assume that after a collision (possibly 
multiple) with a single cube, the molecule leaves the surface and cannot collide with other 
surface cubes. These assumptions hold for a small enough corrugation (sufficiently small 

Q^max) • 

For a molecule falling on a corrugated surface with an incident angle #i nc , the probability 
of hitting the surface depends on the horizontal coordinate X' . This probability is also 
proportional to the cosine of the effective incident angle 8 mc — a: 



Prob(a) oc cos (0 inc — a) = cos 9 inc cos a + sin 6* inc sin a (4) 

1 + tannine tana) . 



cos 6» inc 



vTT tan 2 a 

Substituting Eq. fl3]) and normalizing gives the dependence on X'\ 

Prob(X') = *L l+ ;tanfl inc sin( 27 rX7c ) (5) 
Kc + r 2 sin 2 (2irX'/c) 

where r = 2na/c, and 

K f - * (6) 

J v/l + r 2 sin 2 £ 
o v 
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B. Vibration of surface atoms - simple inclusion of phonons 



The hard cube model [26|, |27|, |29| provides a simple way for including surface atom vibra- 
tion into the collision process. The collision of a molecule with a hard wall of infinite mass is 
replaced by a collision with a hard cube of a finite, but a large mass M moving with velocity 
U. According to the model, the cube oscillates in the direction locally perpendicular to the 
surface. The reason is that for a flat and frictionless cube only the vertical velocity compo- 
nent can transfer energy to the impinging molecule. We slightly modify the original hard 
cube model and let the cube perform harmonic oscillations with a frequency um (the spring 
constant of the oscillator is M(jj 2 M ). The position of the face of the cube, in the unprimed 
coordinate system, is given by: 

Z M {t) = A sin (u M t + (p) , (7) 

and the velocity, therefore, is: 

V z>M (t) = U{t)= u M A cos (u M t + <p) , (8) 

where A and <p are the amplitude and the phase of the oscillator, respectively. These A and 
ip will be chosen randomly, according to the distribution functions, as explained below. 

When treating the collision process, we assume that the mass of the surface atom is much 
larger than the mass of the molecule, i. e. M 3> 2m, so that the collision does not change 
the velocity of the hard cube. Only the velocities of the molecule change, and its total 
energy in the laboratory coordinate frame may increase or decrease. If we assume that the 
oscillating hard cubes are always in thermal equilibrium, then the constant total energy of 
the harmonic oscillator, equal to Mu^A 2 /2, leads to the following probability distribution 
of the amplitudes A: 

where T sur f is the temperature of the solid surface, and ks is the Boltzmann constant. The 
phase (p of the oscillator is uniformly distributed between and 2n. 

The collision condition for a molecule falling on a surface cube can be expressed by 

Z(t)-Z M (t) = K e \ C os6(t)\ , (10) 
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where Z(t) is the time-dependent position of the center of mass of the colliding molecule, 
ZM{t) is given by Eq. (J7]), and 0.5r e cos# is the z-component of the orientation vector r of 
the molecule. 

The translational velocity of the molecule is transformed before the collision according to 
Vz — > Vz — U(t co \), where t co \ is the moment of the collision. This is a usual Galilean trans- 
formation to the coordinate frame moving with the hard cube. Similarly, after the collision 
the translational velocity of the molecule is transformed back by Vz — > Vz + U (t co i) (the 
collision and the change of velocities occur instantaneously in our model). The rotational 
velocity remains unchanged under this Galilean transformation. 



C. The collision of the molecule with a frictionless surface 



The collision process in the local (unprimed) coordinate system can now be treated in 
the same way as in [22]. Here we provide the final results, taken from the Appendix section 



of 22] (after substituting /i = 1 for the homonuclear diatomic molecule.) 

Here and below, the subscripts i and / denote the values of variables before and after 
a single collision, respectively. For a collision with a frictionless cube the parallel velocity 
component is always conserved: Vxj = Vx,%- The colliding molecule can be therefore treated 
as if it had only a Z-component of its center-of-mass velocity. The final perpendicular 



velocity component is: 



T , -cos 2 6^ - 2|v;| sin6> (e v . ■ e v Vi ) 

v Z,f = — — • (11) 

1 + sin 6 



Here e Vi is the unit vector in the direction of the initial rotational velocity Vj, and e Vf - Vi is 
given by: 

e z - cos 9 Tj 
1 sin 6 

where ez is the unit vector in the direction of Z. The final rotational velocity v/ can be 
found using: 

v/ = Vi + \ V ZJ - V z ,i\ sin 9 e Vf _ Vi . (13) 

Before, after and between the collisions (if there are several), the motion of the molecule 
is in free space. The translational motion of the center of mass is: 

Z(t) = Z + V z>0 (t-to) ■ (14) 
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The rotational motion of the molecule is a free rotation with a constant speed vq. The 
plane of rotation is fixed and is perpendicular to the angular momentum vector, which is a 
constant of motion for free rotation. The orientation of the molecule as a function of time 
is: 

r(t) = r cos [u (t - *„)] + ^— sin [u (t - t )} , (15) 

Z Vq 

where vq = |v |, coo is the angular rotational speed of the molecule, and vq = coor e /2. Taking 
the derivative, the velocity is: 

v(i) = -uj t q sin [u (t - t )] + v cos [u (t - t )) . (16) 



III. INTERACTION OF THE MOLECULE WITH AN ULTRASHORT LASER 
PULSE 



Here we briefly summarize the results of the classical model describing the interaction 
of the diatomic rigid molecule with a nonresonant ultrashort laser pulse, in the impulsive 



approximation. A more detailed description may be found in [30 ]. 

The potential energy of the laser pulse interacting with the induced molecular dipole 
moment is given by: 

V(9',<p',t) = -~£ 2 (t) (Aacos 2 ft' + a ± ) , (17) 

where Aa = a» —a± is the difference between the polarizability along the molecular axis and 
the one perpendicular to it, £(t) is the envelope of the electric field of the linearly polarized 
laser pulse, and ft' = ft' (6', (p') is the angle between the molecular axis and the direction of 
polarization of the pulse. We assume that the pulse duration is very short compared to the 
rotational period, so that the pulse can be described in the impulsive (5-kick) approximation. 
We define the dimensionless pulse strength P, as 

Aa 

4h 



/oo 
£ 2 (t)alt . (18) 
-oo 



Planck's constant h appears here and in the remainder of the paper only for the reason 
of unit convention. Under this convention, the quantity HP is the typical value of angular 
momentum transferred by the pulse to the molecule. 

We consider the action of a pulse linearly polarized along some arbitrary unit vector p, 
and determine the vector of the resulting velocity change Av' for a molecule oriented along 
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FIG. 2: (Color online) The suggested experiment, in which the molecular beam falling on the 
corrugated surface at an angle 9 mc = 0° (blue thin arrows) is irradiated by an ultrashort laser 
pulse (thick red arrow), polarized as shown. This laser pulse aligns the molecules and produces a 
time-averaged angular distribution of them, as depicted by the gray cigars. These aligned molecules 
fall on two different slopes of the corrugated surface. They receive a differently directed "kick" on 
the two slopes (small red arrows), and as a result, they scatter to different angles, with a different 
sense of rotation, as shown. 

some direction ro'. The final result is: 

Av' = ^ cos P' (jP- cos #r ') , (19) 

where / = mr 2 e j2 is the moment of inertia of the molecule. 



IV. A SURFACE AS A "TENNIS-RACQUET" SELECTOR OF MOLECULES: A 
CORRELATION BETWEEN THE SCATTERING ANGLE AND THE SENSE OF 
ROTATION 



In our previous publication [22| we suggested a way to produce diatomic molecules, such 



as N 2 , with a preferred sense of rotation using an ultrashort laser pulse and an event of 
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FIG. 3: Probability distributions of the scattering angles of non-rotating (T ro t = 0) nitrogen 
molecules scattered from a frozen corrugated surface (T sur f = 0, a max = 20°), with an incident 
angle of (a) #i nc = 0° and (b) 6 mc = —10°. The corrugation rainbows can be seen at expected 
angles of Rainbow = — #inc ± 2a max . In (b) the left peak is higher than the right peak, because the 
incident molecules have a higher probability to fall on the positive slope of the corrugated surface, 
and to scatter to negative angles. In this Figure, the results of a Monte Carlo simulation of 10 5 
molecules are shown. 

molecular scattering from a flat solid surface [such as Ag(lll)]. We suggested to shine a 
femtosecond laser pulse polarized at +45°/ — 45° to the surface normal on molecules in a 
molecular beam. As a result, the flying molecules will be aligned (on average) along the 
direction of polarization of the pulse. After hitting the flat surface, the molecules will be 
"kicked" by it and receive a preferential rotational velocity of clockwise/counter-clockwise 
sense in the plane containing the polarization vector of the laser pulse and the surface 
normal. 

Our current suggestion utilizes laser-aligned molecules that are scattered from a corru- 
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gated surface, such as LiF(OOl), or Pt(211). The idea is based on the following consideration. 
Consider diatomic molecules that are "shot" towards the corrugated surface with 9\ nc = 0°, 
see Fig. [2J Part of these molecules will fall on a corrugation with negative slope and will be 
scattered mainly to positive scattering angles 6 scaut , while the other part will fall on a corru- 
gation with positive slope and will be scattered mainly to negative scattering angles. Now 
consider irradiating these molecules with an ultrashort laser pulse polarized along the surface 
normal, at 0°, before they hit the surface. The molecules will now be aligned on average at 
0° when hitting the surface. The molecules hitting the negative slope of the corrugation will 
scatter to positive angles and receive a "kick" providing them a counter-clockwise rotation, 
while the molecules hitting the positive slope of the corrugation will scatter to negative angles 
and receive a "kick" providing them a clockwise rotation. Therefore, a correlation between 
the angle of scattering and the direction of rotation should be present in the distribution of 
the scattered molecules. 

When atoms or molecules are scattered from a corrugated surface, rainbows might be seen 



in the scattering angle distribution of the scattered atoms/molecules, see the review 31[ and 



the experiment in [24]]. The angles corresponding to the rainbow scattering are determined 
by the points of inflection of the corrugation function. For a frozen (zero temperature) hard- 
wall surface with a corrugation given by Eq. (El) , the rainbow scattering angles are easily 



shown to be equal to 



3l| 



^rainbow ^inc i 2Q max . (20) 

The simulation results discussed below were obtained using the models of Sections [TT] and 
IIHI and a Monte Carlo averaging. In Fig. [3] we plot the angular distribution of non- rotating 
molecules scattered from a frozen corrugated surface (T sur f = and a max = 20°). In panel 
(a), the incident angle is 6*i nc = 0°, and in panel (b) ^inc = —10°. In (a) the peaks correspond 
to Eq. (]20|) and are symmetric, because the incident molecules hit with equal probability the 
positive and the negative slopes of the surface corrugation. In (b) the peaks also correspond 
to Eq. fT2"0]) and the left peak is higher than the right peak. The reason is that for negative 
incident angles the probability of hitting the positive slope of the corrugation and scattering 
towards negative angles is higher than the same probability for the negative slope. 

In Figs. HI [5] and El we plot the results for a molecular beam of nitrogen molecules falling 
on the corrugated surface (a max = 20°) at several incident angles and at a constant speed 
of Vq = 300 m/s. The molecules in the beam have a rotational temperature of T rot = 5K, 
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FIG. 4: (Color online) A two-dimensional probability distribution of the scattered molecules. The 
horizontal axis represents the scattering angle, and the vertical axis is for the angular momentum 
component along the Y'-direction, in units of H. Hot (red) colors represent high probability, while 
cold (blue) colors represent low probability. Before the scattering, the molecules were "kicked" by 
a laser pulse polarized at 0°, along the Z'-axis, with a pulse strength of P = 10. The angle of 
incidence is 9 mc = 0°. A clear correlation between the scattering angle and the sense of rotation 
is seen. In (b), a distribution of the scattering angles of the molecules is shown. The incident 
nitrogen molecules {m = 14 a. u.) have a speed of Vq = 300 m/s and a rotational temperature of 
Trot = 5K. The corrugated surface is composed of platinum atoms (M = 195 a. u.) vibrating at a 
frequency of ojm = 10 13 Hz, while the amplitudes of vibration are distributed according to Eq. Q 
with T sur f = 300 K. The corrugation strength of the surface is a max = 20°. In this Figure and in 
Figs. [5] and [6] the results of Monte Carlo simulations of 2 x 10 5 molecules are shown. 
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FIG. 5: (Color online) The same distribution as in Fig. d] is plotted, except the incident angle that 
is Bine = —10° here. The slight rotation of the molecular beam (or the surface) leads to a change 
of the relative intensity of the two rainbow peaks. 

typical for molecular beam experiments. The surface is composed of cubes of mass M = 
195 a.u., corresponding to platinum atoms, with oscillation frequency of um = 10 13 Hz. The 
amplitudes of the vibrating cubes are distributed according to Eq. (j9]), with a temperature 
of T surf = 300 K. 

In Fig. HI we "shoot" the molecules towards the corrugated surface at normal incidence 
(#inc = 0°), "kick" them with an ultrashort laser pulse polarized along the same vertical 
direction, and let them scatter from the corrugated surface In Fig. H^a) we plot a two- 
dimensional histogram, where the vertical axis represents the angular momentum component 
Lyi perpendicular to the plane of the expected molecular rotation, and the horizontal axis 
represents the scattering angle 6 sca ,t- A clear correlation can be seen between the angle of 
scattering and the sign of the angular momentum component. Molecules with positive Ly> 
are scattered to negative angles, as expected. 
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The distribution of scattering angles [which is also shown in the two-dimensional plot 
in Fig. H(b)] exhibits a double rainbow form, as in Fig. [3j However, the location of the 
peaks is at smaller values of |# scat | than the values given by Eq. (120]) . This happens for 
two reasons. First, the surface has a non-zero temperature, and the vibrating surface cubes 
transfer translational velocity to the molecules in the Z-direction (the corrugation is rather 
small and all the cubes vibrate mainly along the Z direction). This transfer of velocity moves 
the scattering angle closer to the surface normal. Second, the laser pulse provides a high 
rotational energy to the molecules. During the collision this high rotational energy is partly 
transferred to a translational energy. Because the Vx velocity component is conserved, most 
of the energy goes to the Vz component, which, because the surface corrugation is small, 
again, moves the scattering angle closer to the surface normal. 

The same kind of two-dimensional distribution for 8- mc = —10° is shown in Fig. |5j Here 
the two rainbow peaks move to different angles, as seen qualitatively from Eq. ( )20|) . and their 
relative intensity changes, as explained above. This suggests a way to control the intensity 
of molecular beams that are scattered from a solid surface with a preferred sense of rotation. 

Going back to the case of 9- mc = 0° and changing the polarization of the laser to 90°, we 
obtain Fig. [6j The correlation between Ly and (9 scat is now inverted, as could be expected. In 
addition, the effect now is much weaker than for the case of 0° polarization of the laser pulse, 
in Fig. HI The reason is that now the molecules that are falling on the surface are aligned 
parallel to the surface. Because the surface is almost flat (a max = 20°), these molecules can 
collide twice with the surface, when the second collision almost cancels out the rotational 
excitation caused by the first collision. 



V. A MODEL OF MOLECULAR SCATTERING FROM A SURFACE WITH 
FRICTION 



In this Section we develop a model and analyze the scattering of a molecule from a surface 
that is "corrugated" at a length scale much smaller than the size of the molecule. This 
"corrugation" is treated as a friction force acting parallel to the flat surface. The prototype 
for our model is the model that was used in order to describe qualitatively some effects in 
the scattering of nitrogen molecules from the Ag(lll) surface 25]. This silver surface is 
practically flat, however, the measured effects in [25| suggest that in-plane tangential forces 
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FIG. 6: (Color online) The same distribution as in Fig. H] is plotted, except the direction of the 
polarization of the laser pulse, that is 90° here, or along the X'-axis. Now, the incident molecules fall 
almost parallel to the weakly corrugated surface. This leads to double-collisions of the molecules, 
where the second collision almost cancels out the rotational excitation of the first collision. This 
still leaves some correlation, which is, however, much weaker than in Fig. 01 

are present during the collision process. These in plane tangential forces can be created 
by different physical mechanisms, such as: tangentially directed phonons, electron-hole pair 
creation, etc. 



Similar to Sec. [TTJ we extend our model introduced in reference 22]. In this Section, we 
add a friction force between the colliding atom and the flat surface, in a similar way to 25 ]. 
The friction force is exerted parallel to the surface, its magnitude is proportional to the local 
velocity of the colliding atom and its direction is opposite to this local velocity direction. In 
subsection IV B\ the model is developed for a frozen surface, and for a two-dimensional planar 
rotation of the molecule. In subsection IV C\ the model is extended to a three-dimensional 
rotation. The extension to a finite-temperature surface is done as in subsection III Bt 
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A. Definition of the coordinates 

The model we use in this and the following Sections treats the surface as a flat surface 
with friction, rather than a corrugated wall (later in subsection I VI C[ the surface is composed 
of vibrating hard flat cubes). The variables and the parameters are similar to those of Sec. 
Hlj with few exceptions. The (X, Y, Z) and the (X', Y' , Z') coordinate systems are now 
identical one to another. The coordinate system from this point on is (X,Y,Z), as shown 
in Fig. [7J Because the frictional forces can act also in the Indirection, the center-of-mass 
velocity should include now three components: V = (Vx, Vy, Vz\ The molecular orientation 
is r, and the rotational velocity is v. The mass of the molecule is 2m. The incident angle 
is #mc = arctan (Vx,o/Vz,o) (initially Vy is always chosen as zero). The point of collision is 
denoted by O. The colliding atom is denoted by A, and the other atom is denoted by B. 
It is assumed that the special case when both atoms collide simultaneously with the surface 
never occurs. 

B. A two-dimensional model of molecular scattering from a surface with friction 

If the translational and rotational motion of the molecule is restricted to the XZ-plane, 
the Y coordinate and the Vy velocity component are always equal to zero. The orientation 
of the molecule in the XZ-plane can now be defined by a single angle 9 that increases in a 
clockwise direction from the direction of the Z-axis and ranges from to 27r. The rotational 
velocity is now defined by a single variable v = |v^| = |vb|, equal to ur e /2, where u is the 
angular velocity of the molecule. 

1 . Friction force 

At the moment of the collision, two forces act on the colliding atom A: a normal force 
F n and a friction force Ff r , see Fig. [7J We assume that F n is a conservative force, that its 
functional form is not known. However, we assume that Ff r is a nonconservative force that 
has a known functional form, given below. We also assume that the time interval during 
which the forces are applied is very short compared to the time it takes to the variables X, 
Z and 9 to change significantly, such that the forces are practically instantaneous. 
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FIG. 7: (Color online) Diatomic molecule (a dumbbell) colliding with a flat, rough surface. The 
molecule has translational velocity V, with X and Z components Vx and Vz, respectively. The 
rotational velocity of atoms A and B is vj± and vb, respectively, where |v^| = |vg| = v. The 
position of the center of mass is given by the coordinates X, Y and Z, and the orientation vector 
of the molecule is r. The forces acting on the atom A at the point of contact O are shown in red: 
F n is the normal force, and Ff r is the friction force, given by Eq. (|2ip for the two dimensional case, 
and by Eq. (|3T[) for the three-dimensional case. 

As explained above, the friction force is given by: 



where 7 is the friction coefficient, ex is a unit vector in the X direction, and 5(t) is the 
Dirac delta function in the time variable. The time t = denotes the moment the molecule 
arrives to its turning point during the collision with the surface. It is assumed here and 
below, that 7 -C 1. The meaning of this assumption will be clear in the next paragraph. 

2. Solution of equations of motion and conservation of angular momentum 

Momentum change in the X-direction. In the following derivations, we assume that 




(21) 
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the orientation 9 of the molecule is in the range between and 7r/2. Later, we extend the 
final results for any angle. Newton's equation of motion for the center of mass motion in 
the X-direction is: 

= -2m 7 [V x {t) - v{t) cos0(t)] 8{t) . (22) 

(Jj V 

After integrating over the time of the collision, we obtain: 

Vx ~ Vx = -7 [Vx(0) - v(0) cosO] , (23) 

where the superscripts — and + denote the values of the variables immediately before and 
after the collision, respectively. We assume that during the collision the position and the 
orientation of the molecule do not change, while the velocities do change. We assume that 
the velocities at the turning point Vx(0) and t>(0) are equal to the velocities before the 
collision V% and v~, respectively. This assumption is valid, because our basic assumption is 
that 7 1, and this means that the change in Vx as the result of the friction force is small, 
compared to Vx itself. We obtain: 

V£ = Vx-7(Vj-v-\<x*0\) , (24) 

where the absolute value makes the expression valid also for values of 9 higher than 7r/2. 
Defining the effective speed of the colliding atom before the collision as: 

Vx& = Vx-\case\v- , (25) 

Eq. (|24|) can be written as: 

n = V x - 7V£> • (26) 

Angular momentum conservation with respect to the origin O. Taking the point 
O to be the origin (see Fig. [3), it can be seen that the lever arms of the forces F n and Ff r are 
zero. This leads to the conclusion that the collision process preserves the angular momentum 
of the molecule, measured with respect to the origin O. In addition, the angular momentum 
of atom A is zero, so the conservation law is: 

L + B = L B , (27) 

or, explicitly: 

cos9V^ - sin9V^ + v + = cos9V x -sin9V z +v~ , (28) 



20 



similar to Eq. (4) in 22]. 



Energy change. The total molecular energy is the sum of the translational energy of 
the center of mass and the rotational energy with respect to the center of mass, and is given 
by: 

£=i(2m)(yi + V|) + i(2my . (29) 

The rate of change of the energy during the collision, or the power, is given by: 

dE 

— = (F n + F fr )-(V + v A ) . (30) 

Integrating over the collision time, noticing that F n is a conservative force that does not 
contribute to the change of energy and using Eq. (12~T1) . we obtain: 

E + -E- = -2m 7 [V x (0) - v(0) cos0] 2 . (31) 



Using Eq. (|29|) . assuming that the velocities at the turning point are equal to the velocities 
before the collision and performing some rearrangements, we arrive at: 

(V^) 2 + (^ + ) 2 + K) 2 = (l-27)^^ (32) 
3. The final velocities 

The velocity V+ is given by Eq. ((26]). Using Eqs. fl26]), (J2HJ) and $2$ the values of V+ and 

v + can be found using tedious but straightforward algebra. We solve a quadratic equation, 

and the sign before the square root is chosen such that the following equations for 7 = 

are consistent with Eqs. ( fill and ( {TBI . At the end, we generalize the final results for the 

orientation angle 9 in the range from to 2ir, by adding an absolute value [as in Eq. fl2~4"j) ]. 

and a correct sign, where necessary. We also neglect terms proportional to 7 2 , while leaving 

terms proportional to 7, for consistency with the assumption 7CI. The final result is: 

cos 2 9 V% + 2 sgn(tan 9) | sin 9\v~ + 7 sin 29 V^ cS 
Vz = 1 + sin 2 9 '~ • (33) 

Using Eqs. (l26j) and ( J28l) . the result for v + , expressed using V£ is: 

v + = + sgn(tan^)] sin^|(U+ - V z ) + 7 | cos9\V XcS , (34) 



or, after substituting Eq. (|33|) : 

cos 2 9v~ - 2 sgn(t&n 9)\ sin 9\Vz +7] cos 3 ^] V^ eff 
V ~ 1 + sin 2 9 
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(35) 



In the equations above, sgn(a) denotes the signum function, defined as sgn(a) = a/\a\. 

The free motion of the molecule before, after, and in between the collisions is according 
to Eqs. (JHJ), (EES]), (TToT) . Equations (TTo"|) and (TT6"j) for a two-dimensional rotation reduce to: 



The condition for collision of the molecule with the surface is given by Eq. ( ITOj) . where 
Zhiif) = for a frozen surface. 

C. A three-dimensional model of molecular scattering from a surface with friction 

Now we extend the model to the case of a three-dimensional rotation. By our convention 
(see Fig. [7]), at the moment of collision the atom touching the surface is denoted by A, while 
the other atom is denoted by B. It follows that at the moment of collision the orientation 
of atom A is given by r A = — r, if the ^-component of the orientation vector r is positive, 
and it is given by ta = r, if the ^-component of the orientation vector r is negative. 

The rotational velocity of the molecule is now given by the vector v, which is changing 
with time, as opposed to the scalar quantity v = |v|, which is constant for a freely rotating 
molecule. The components of the velocity v are given by (v x , v y , v z ), where the (x,y,z) 
coordinate system is fixed to the center of mass of the molecule and is parallel to the 
(X, Y, Z) system. 

In the three-dimensional rotation case, the velocity component Vy, that is initially zero, 
can become nonzero, because of the friction forces acting in the ^-direction. 

1. Friction force 

The friction force is now proportional to the local velocity component of atom A parallel 
to the surface, and is given by: 




(36) 



F fr = -2m 7 [e x (V + v A ) ■ e x + e Y (V + v A ) • e Y ] S(t) . 



(37) 
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2. Solution of equations of motion and conservation of angular momentum 

Similar to Eq. 021]) we obtain for the velocity components V x and Vy after the collision: 

n = V X -1 {V x + v~ Ax ) , (38) 

and 

V+ = Vy - 7 (Vy + v- A>y ) . (39) 
Finding the energy change according to Eq. (I3"U|) . we obtain, similar to Eq. (IBTi) : 

E+ - E- = -2m 7 [(V x + v~ AjX f + (Vy + v^f] . (40) 

In the last equation E denotes the total molecular energy, equal to: 

E = ~(2m) (Vl + V* + Vl + v% x + v\ y + v% z ) . (41) 

The equation for angular momentum conservation, similar to Eq. fl28|) . is now expressed 
in vector form: 

r B x (V+ - V- + v+ - vi) = , (42) 

or, changing the orientation and rotational velocity vectors from the non-colliding atom B 
to the colliding atom A: 

r A x (V- - V+ + v+ - v^) = . (43) 

3. The final velocities 

Solving together Eqs. (155]). fl39]) . f]40]) and ff43]) (Mathematica 7 software was used), and 
neglecting terms proportional to 7 2 or a higher power, we obtain the following expression 
for V£: 

sin 2 9(V Z —v A ~ z ) + 0.5 sin 26 (v± - jV±) - sgn(iT) (VX - 2~f sin 9 V^ R~ / VX) 

y+ — ! ^ L 

z ~ 1 + sin 2 ^ 

(44) 

where 

v± — v A,x cos f + v A,y si 11 V 9 5 (45) 
V± = (V x + v A>x ) cos <p + (V Y + v AtV ) sin Lp , (46) 
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and 



R = sm9v x - -cos6{V z -v A>z ) , (47) 



A = 2R 2 + ^(l + sm 2 6)(V z + v Az ) 2 . (48) 



Eq. ( 144")) reduces to Eq. ( 133|) in the special case of rotation in the XZ plane. In addition, 
we checked numerically that Eq. (jUJ) reduces to Eq. f JTT]) for 7 = 0. From a vector diagram 
drawn according to Eq. ( 143|) . and from the fact that is always perpendicular to r A , we 
find an expression for vj[, similar to Eq. f Tl~3|) : 

v a = v a + I v+ " v I sin a e v +_ v - , (49) 

where 



and 



sin a = y 1 - (r A • e) , (50) 
e- (r A -e)r A 

e v +_ v - = : , (51) 

A A sma 



V" - V + 

(52) 



|V- - V+| 

The last equations reduce to equations from Sec. Ill Cl if e is replaced by ez- This can be done 
in the case of scattering from a frictionless surface, because in such a collision the Vx and 
Vy components are conserved, and the collision can be treated as a collision perpendicular 
to the surface, leading to e = e z . 



VI. SCATTERING OF UNIDIRECTIONAL ROTATING MOLECULES FROM A 
SURFACE WITH FRICTION 

In the current Section, we investigate the influence of the direction of rotation of molecules 
on their scattering from a surface with friction. In subsections IVI Al and IVI Bl we consider 
some idealized cases of non-rotating and rotating molecules randomly oriented in the XZ- 
plane, that are scattered from a frozen surface with friction. After clarifying the "clean" 
effects seen for two-dimensional, zero temperature surface conditions, we treat the three- 
dimensional molecular rotation, finite temperature vibrating surface case, with addition of 
laser pulses, in subsection IVI CI 
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FIG. 8: (Color online) Rotational velocity of the scattered molecules is plotted as a function of the 
initial orientation angle of the molecules in the XZ-pl&ne. Initially, the molecules are non-rotating, 
and they impinge on a frozen surface with an incident angle of 6 mc = —45°. The solid (blue) line 
represents scattering from a frictionless surface (7 = 0), while the dashed (red) line is for scattering 
from a surface with friction (7 = 0.1). 

In all the following subsections, the molecules collide with the surface with friction at 
an incident angle of —45°. In subsections IVI Al and IVI Bl the values of the velocities are all 
given in units of Vb, such that Vx,o = \Vz,o\ = l/\/2. 

A. Scattering of non-rotating molecules — a two-dimensional case 

In order to check the two-dimensional model of Sec. IV B\ we inspect the final scattering 
angles and the final rotational velocities for initially non-rotating molecules, oriented in the 
XZ-plane, that scatter from a frozen surface with friction. In Fig. [8] we plot the dependence 
of the final rotational velocity of the scattered molecule on the initial orientation angle 
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FIG. 9: The distribution of the scattering angles of 10 5 molecules scattered from a frozen surface 
(^surf = 0)- The molecules are non-rotating and are randomly oriented in the XZ-plane before the 
scattering, and they impinge on the surface with 9 inc = —45°. (a) The surface is frictionless with 
7 = 0, and the scattering is mainly specular, similar to scattering of atoms, (b) The surface is rough 
with a friction coefficient 7 = 0.1, and the scattering is mainly towards angles lower than specular, 
because the friction of the surface slows down the motion of the molecules in the X-direction. 

$o in the XZ-plane. The solid (blue) curve represents the case of a frictionless surface 
(7 = 0), while the dashed (red) curve is for a surface with friction (7 = 0.1). Three regions 
can be identified in the curves of Fig. [SJ the first region of positive rotational velocities 
(0° < 6*o < 45°), the second region of negative rotational velocities (135° < 9 < 180°), and 
the third region of the sharp transition between the first two regions, with relatively small 
rotational velocities (45° < #0 J$ 135°). The first two regions represent a scattering process 
by a single collision of the molecule with the surface. It is seen that in the first region 
the friction speeds up the rotation caused by the collision, while in the second region the 
friction slows down the rotation. This asymmetry of the friction force can be understood 
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FIG. 10: The distribution of the scattering angles of 10 5 molecules scattered from a frozen surface 
C?surf = 0)- The molecules are rotating clockwise in the XZ-plane before scattering, with vq/Vq = 2, 
and they impinge on the surface with 9{ nc = —45°. (a) The surface is frictionless (7 = 0), and the 
molecules scatter mainly to subspecular angles because the energy is transferred from rotational 
to translational degrees of freedom, during the collision, (b) The surface is rough with a friction 
coefficient of 7 = 0.1. The scattering peak moves towards a slightly higher angle than in (a) because 
of the frictional force, as explained in the text, and the width of the distribution is significantly 
higher than in (a). 

if the strong normal force F n and the weak friction force Ff r are analyzed for the molecule 
hitting the surface, see Fig. |7J The normal force gives some velocity to the molecule, which 
is symmetric but opposite in sign for molecules inclined by the same angle to the right 
(0° <0 < 45°) and to the left (135° < 6 < 180°). This can be also seen from the second 
term on the right of Eq. ( I34j) . On the other hand, the friction force acts always in the same 
direction on the molecules falling on the surface from the left (negative incident angles). 
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This way, the friction speeds up the rotation, caused by the normal force, of a molecule 
inclined to the right, however, it slows down the rotation of a molecule inclined to the left. 
This can be seen from the third term on the right of Eq. (1341) , which has a positive sign for 
all the angles 9, for non-rotating molecules. 

The third region of the curves in Fig. [S] (45° < 9q < 135°) is a region where the scattering 
process contains double collisions. Here, the rotational velocity given to the molecule by the 
first collision is approximately canceled out by the second collision. 

In Figure IHa) we plot the scattering angle distribution of non-rotating molecules in 
the XZ -plane scattered from a frictionless surface. The majority of the molecules scatter 
towards the specular angle around 45°, as expected from Fig. [SJ However, from Fig. Mjo) 
we see that when friction is added, two peaks appear in the distribution, and they move 
to smaller scattering angles. The shift to smaller angles can be understood as an effect 
of the rough surface that "slows down" the translational motion of the molecules in the 
X-direction. The peak at 42° in Fig. Mjo) corresponds to molecules that collide once with 
the surface, while oriented at and near 6 = 0°. The higher peak at 39° corresponds to 
molecules oriented around 9 = 90° and colliding twice with the surface. The molecules 
colliding twice lose more of Vx than the molecules colliding once, and scatter, as a result, 
to smaller scattering angles. The location of the peaks for non-rotating molecules can be 
found according to: 

# S cat, P eak = arctan — — = arctan — — , (53) 

Vzj \ V Zfl\ 

where n is the number of collisions, that can be equal to 1 or 2. 



B. Scattering of unidirectional rotating molecules — an idealized two-dimensional 

case 

Here we examine the ideal scenario of molecules in the XZ-plane, all rotating unidi- 
rectionally at a constant angular speed, and colliding with a frozen surface with friction. 
In Sec. IVI CI we extend this scenario to more realistic conditions: molecules set in a three- 
dimensional rotational motion by femtosecond laser pulses, colliding with a vibrating surface 
at a finite temperature. 

Unidirectionally rotating molecules colliding with a surface with friction can be viewed 
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FIG. 11: The distribution of the scattering angles of 10 5 molecules scattered from a frozen surface 
(^surf = 0). The molecules are rotating counter-clockwise in the XZ-pl&ne before scattering, with 
v o/Vo = — 2. They impinge on the surface with 9{ nc = —45°. (a) The surface is frictionless (7 = 0). 
This Figure is the same as Fig. [TUfa). and is given here for comparison with panel (b). (b) The 
surface is rough with a friction coefficient of 7 = 0.1. The scattering peak moves towards a lower 
angle than in (a) because of the frictional force, as explained in the text. 

as a rotating wheel falling on the ground. Inverting the direction of rotation of the wheel 
should change the scattering process, as explained below. 

Imagine a wheel falling on the rough surface from the left, while rotating clockwise. Let 
us suppose that the rotational velocity of the wheel is high enough, such that the effective 
velocity of the wheel at the point of contact is directed to the left. In this case, the friction 
force that is exerted by the surface at impact is directed to the right. This friction force will 
cause the wheel to scatter towards scattering angles higher than at the absence of friction. 
On the other hand, for a wheel that rotates counter-clockwise, the friction force created at 
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the point of contact is directed to the left. This friction force will cause the wheel to scatter 
towards scattering angles lower than at the absence of friction and closer to the surface 
normal. 

Now consider molecules initially rotating in the XZ-plane with a constant speed of 
v o/^o = +2 for the clockwise rotation, and of Vq/Vq = —2 for the counter-clockwise ro- 
tation, and colliding with the surface with friction. The effect of the shift in the scattering 
angles is visible if we compare Figure ITuT b) for the clockwise rotation, to Figure [TTT b) for 
the counter-clockwise rotation. For comparison, the corresponding scattering angle distri- 
butions for rotating molecules scattered from a frictionless surface are shown in Figs. [TUT a) 
and [TIT a). These two Figures are identical, because the sense of rotation does not change 
the scattering process for a frictionless surface. A clear deviation of the scattering angle in 
opposite directions for clockwise and counter-clockwise rotation is seen in Figs. [TOT b) and 
[TTT b). However, this deviation is lower for the clockwise rotating molecules, than for the 
counter-clockwise rotating molecules. The reason for this is the difference between the two 
cases in the local velocity of the colliding atom. The local velocity parallel to the surface 
of the colliding atom is given by Vx,o + va,x- Let us consider a rotating molecule oriented 
at 0° colliding with the surface. In the clockwise rotation case the local velocity of this 
molecule is Vx,o ~\ v o\ — —1-3, while in the counter-clockwise rotation case the local velocity 
is Vx,o + l^o I = 2.7. The higher (in absolute value) local velocity in the counter-clockwise 
rotation case leads to a higher friction force, that results in a higher change of the scattering 
angle. 

C. Scattering of unidirectional rotating molecules — a realistic three-dimensional 
case, with laser pulses 

After demonstrating the difference between the scattering of molecules rotating in an 
opposite sense for the idealized two-dimensional case, we turn to more realistic conditions, 
involving the use of laser pulses. In this subsection, we model a molecular beam of nitrogen 
molecules (atomic mass m = 14 a. u.) falling on the surface with friction at an angle of 
incidence of —45° and at a constant speed of Vq = 300 m/s. The molecules in the beam 
have rotational temperature of T rot = 5 K, which is typical for molecular beam experiments. 
The surface is composed of cubes with M = 108 a. u., corresponding to silver atoms, with 
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oscillation frequency of % = 10 13 Hz, corresponding to the Ag(lll) surface [32|. See Sec. 
Ill Bl for the definition of the parameters. The surface is at room temperature, T sur f = 300 K, 
and has a friction coefficient of 7 = 0.1. 



Before hitting the surface, the mo 
following the scheme introduced in 



ecules are illuminated by two consecutive laser pulses, 



331 ] . and experimentally realized in |34|. See Section 



II III for the rotational velocity change caused by a single laser pulse. The first linearly 
polarized ultrashort laser pulse causes transient alignment of the molecules. The second 
pulse, hitting the molecules at the moment of maximal alignment, and polarized at 45° 
to the first pulse, causes the molecules to rotate preferentially in the plane defined by the 
polarization directions of the two pulses. The sense of rotation (clockwise vs. counter- 
clockwise) is defined by the relative orientation of the polarization directions of the first and 
the second pulse. We want to create a rotation in the XZ-plane, and choose the polarization 
direction of the first pulse to be along the Z-axis. This fixes the polarization direction of 
the second pulse to ±45° from the Z-axis, where the plus sign corresponds to the ignition of 
the clockwise rotation in the X Z-plane, while the minus sign ignites the counter-clockwise 
rotation. The pulse strength chosen is Pi = 5 for the first pulse and P2 = 10 for the second 
pulse. The value of P2 is higher than that of Pi, unlike in [33], in order to create higher 
rotational velocities. 

In Fig. [T2]we show two distributions of scattering angles, corresponding to the simulations 
with the parameters above. The dashed line (blue) distribution corresponds to molecules 
set in a clockwise rotation by the laser pulses, while the solid line (black) distribution corre- 
sponds to the counter-clockwise rotating molecules. A distinction of about 5° can be made 
between the peaks of the two distributions. The distinction between the oppositely rotating 
molecules is not as sharp as in the idealized case of Figs. [10] and [TTJ and the two distribu- 
tions in Fig. [12] have a significant overlap. However, at the left edge of the distributions a 
significant contrast exists between the numbers of oppositely rotating scattered molecules. 
As an example, between the angles of 10° and 15° to the surface normal, the probability to 
find counter-clockwise rotating molecules is about one order of magnitude higher than the 
probability to find clockwise rotating molecules. 
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FIG. 12: (Color online) Two distributions of the scattering angles of molecules set in unidirectional 
rotation by two laser pulses, as explained in the text. The dashed (blue) line corresponds to 
clockwise rotating molecules, while the solid (black) line corresponds to counter-clockwise rotating 
molecules. In this Figure the results of a Monte Carlo simulation of 10 5 molecules are shown. 

VII. CONCLUSIONS 

We introduced two simple classical models, that are both based on the hard cube model, 
and both treat the diatomic molecule as a dumbbell. The models were used to describe 
the laser-induced modification of the molecule-surface scattering process. The first model 
treated the surface as a frictionless surface with a corrugation geometry. The second model 
treated the surface as flat, with tangential friction forces present. 

Using these two models, two potential schemes of laser-controlled surface scattering were 
introduced and investigated. In the first scheme, the molecules were aligned parallel or 
perpendicular to the surface by a laser pulse, before hitting the surface. These molecules 
scattered mainly in two different directions, according to the surface corrugation, and the 
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angle of incidence. The molecules in these two directions were found to be rotating in an 
opposite sense. This difference in angular momentum between the molecules can, in prin- 
ciple, be detected using resonance-enhanced multiphoton ionization (REMPI) spectroscopy 



25|, 



34| . In the second scheme, the molecules were prepared in a state with a specific sense 



33], can provide new means for 



of rotation before the scattering from the surface. It was shown that because of the friction 
with the surface the scattering angles of the rotating molecules depend on their initial direc- 
tion of rotation. This result, combined with suggestions in 
separation of molecular mixtures of different species (such as molecular isotopes, or nuclear- 
spin isomers). It is our hope, that this work will arouse interest in the femtosecond lasers 
community, as well as in the molecule-surface scattering field of research. 
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